Structure and properties of forsterite-MgSiO3 liquid interface: molecular dynamics study

The mechanical properties of partially molten rock, such as their permeabilities and viscosities, are important properties in geological processes. We performed molecular dynamics simulations in terms of structures and diffusivities in forsterite-MgSiO3 liquid interfaces to obtain the nanoscale dynamic properties and structure of the interface. The characteristic structure of the forsterite-MgSiO3 liquid interfaces was observed in the simulations. In the layered structure of the altered surfaces, Si-rich and Mg-rich layers exist alternately in the vicinity of the crystal-liquid interfaces. The layered structure might be formed by the strength difference between Si-O covalent bonds and Mg-O ionic bonds. The difference in the layered structure, indicated by the thickness of the MgSiO3 liquid film, might be caused by the difference in degrees of freedom of the configuration in the liquid film. The two-dimensional diffusivity of oxygen atoms parallel to the interface is controlled by two factors. One factor is the thickness of the liquid film, which decreases oxygen diffusivity with decreasing film thickness. The other is the composition of the sliced layer, where oxygen diffusivity increases with increasing Mg/Si ratio. The effect of the crystal-liquid interface found in this study is negligible in texturally equilibrated rocks. However, the interface can affect the melt flow in deformed samples because a grain boundary melt film with a thickness of several nanometers exists stably in deformed partially molten rock.